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Structural properties of a finite number (TV = 2 — 20) of point charges (classical electrons) confined 
laterally in a two-dimensional two-minima potential are calculated as a function of the distance (d) 
between the minima. The particles are confined by identical parabolic potentials and repel each 
other through a Coulomb potential. Both ground state and metastable electron configurations are 
discussed. At zero distance previous results of other calculations and experiments are reproduced. 
Discontinuous transitions from one configuration to another as a function of d are observed for 
7V = 6,8, ff,f6, 17, 18, 19. 

PACS numbers: 73.22.-f,36.90.-f f,61.46.+w 
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I. INTRODUCTION 

Quantum dots (sometimes called artificial atoms) are 
nanoscale semiconductor structures where a small nuE 
ber of electrons are confined into a small spatial region.l 
The electron motion is usually further restricted to two 
dimensions. There is strong theoretical evidence for the 
existence of a limit where the electron system crystallises 
to Wigner molecules, which is seen as the localisation of 
the electron density_ajpjuiadjpositions that minimise the 
Coulomb repulsion.tlQdu'D'ErEl In the limit of weak con- 
finement (low density) or a very strong magnetic field the 
quantum effects are quenched or obscured and the classi- 
cal electron correlations start to dominate the properties 
of the system. The ultimate limit is a purely classical 
system where only the Coulomb repulsion between the 
electrons defines the ground state. The problem reduces 
to finding the classical positions of electrons (which de- 
pend on the forms of the confining and the interaction 
potentials) that minimise the total eneroz:.! 

There is growiHg-iiiterest in calculati: 
and measuringllZrll3il3 the properties of coupled quan- 
tum dots. Due to the 2D nature of quantum dots the 
two-atom system is different whether the quantum dots 
are coupled in the plane in which the electrons are con- 
fined [laterally coupled) or in the perpendicular direc- 
tion {vertically coupled). Especially for laterally cou- 
pled quantum |-dpts—jj«Jy a limited number of studies 
have appeared.OllSllalla Classical studies serve as a good 
starting point for more demanding quantum mechanical 
calculations. Moreover, the study of classical electrons 
in vertically coupled artificial atoms has revealed inter- 
esting structural transitions in the ground state electron 
configipations as a function of the distance between the 
atoms .c3 

Apart from quantum dots in the classical limit the 
point charges in 2D can be used to model also other 
physical systems. Examples include vortex lines in su- 
perconductors and superfluids and electrons on the sur- 
face of liquid He (see Ref. ^ and references therein). 
In the theoretical field, the ground state configurations 
of a confined classical 2D electron system have been 
studied in the case of a single artificial atom in Refs. 
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61)2 7| and for the vertically coupled ar- 
molecule as a function of the inter-atom 
Recently, also some experimental 
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tificial atom 
distance in Ref. 
studies of 2D confined charged classical particle systems 
haV|e-Ap|peared to reflect the classical cluster patterns in 
2DMEr 

Classical point charges in a two-dimensional infinite 
plane crystallise into a hexagonal lattice at low temper- 
atures. Parabolic confinement in the artificial atom, on 
the other hand, favours circularly symmetric solutions. 
The ground state configuration is thus determined by 
two competing effects, circular symmetry and hexagonal 
coordination, thus resulting in non-trivial particle con- 
figurations. The reported configurations of the electron 
clusters in a single artificial atom do not all agree be- 
tween different studies. The differences can be partly 
explained by the different forms of confinement and in- 
teraction potentials. However, when the number of par- 
ticles, N , confined in the atom is one of the following 
N = 2 — 5,7, 10, 12, 14, 19 all results are in agreement 
while differences appear for = 6, 8, 9, 11, 13, 15 — 18, 20 
(for TV < 20). 

In this paper we consider two laterally coupled artifi- 
cial atoms and classical electrons in the molecule. The 
changes in the ground state electron configurations are 
studied for N — 2 — 20 electrons in the molecule as 
the inter-atom distance is changed. The energies of the 
metastable states are also calculated at different dis- 
tances and their role in the structural transitions in the 
ground state electron configurations is discussed. We also 
reproduce electron configurations of the single parabolic 
artificial atom. The differences between different calcula- 
tions and experimental results are discussed in the limit 
of single atom. 



II. MONTE CARLO SIMULATION 



The classical electrons in the artificial atom molecule 
are modelled with the Hamiltonian 
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H = -m*ujl ^min[(r- -d/2)2,(r, + d/2)2] 
1=1 
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Each of the N electrons is described with coordinates 
•fi = {xi,yi) in two-dimensional space. The harmonic 
confinements are positioned symmetrically around the 
origin with distance d between the minima, m* is the 
electron effective mass, the confinement strength and 
e the dielectric constant. We measure the energy in meV 
and distance in A. The confinement strength was set to 
Hloq = 3 meV and typical GaAs parameters were cho- 
sen to the effective mass and the dielectric constant: 
m* = 0.067 rrie and e — 13. The calculated energy 
values and distances could be scaled to correspond to 
different values of ujq, to*, e, but changing the parame- 
ters also changes the effective distance between the min- 
ima, d, and then the minimum energy configuration may 
not be the same anymore. Therefore we have have one 
significant parameter in the system, d, which scales as 
oc (TO*cj2e)-i/3. 

The minimum energy as a function of the positions of 
the particles, Etot — min i?(ri, ttv), is solved with 
a standard Metropolis Monte Carlo methodic starting 
from a randomly chosen initial electron configuration 
ri, Tat. The accuracy and simulation time needed with 
the Metropolis algorithm was found to be well sufficient 
for the current problem. We compared the calculated en- 
ergies in the limit of a single artificial atom to those given 
in Ref. and the results were found to be in complete 
agreement within the given accuracy. 

In the simulations we choose four different distances 
between the atoms and perform 300 test runs at each par- 
ticle number (A^ = 2 - 20) and distance {d = 0, 200, 600 
and 1000 A). In addition to minimum energy configura- 
tions we also obtain metastable states that are higher in 
energy compared to the ground state. 

When the ground and metastable states are obtained 
at d ^ 0,200,600, 1000 A we study the structural tran- 
sitions between ground state electron configurations at 
the intermediate distances. The electron configurations 
obtained from the fixed d calculations are taken as an 
input to Monte Carlo minimisations where the attempt 
step is set so small that the electron configuration can- 
not change to another. Then the distance is slightly al- 
tered (d — > d ± 1 A) and a new energy with slightly 
modified positions is calculated for the configuration de- 
fined by the input. The calculated new configuration is 
taken as an input to the next calculation with a new dis- 
tance between the atoms, and so all distances between 
d = 0, 200, 600, 1000 A are well sampled. However, it 
may happen that a configuration becomes unstable as 
the distance is changed. In that case the simulation con- 
verges to some other stable configuration, which can be 



seen as a sudden jump to a new energy value in the E{d)- 
plots. The energies of all states are studied as a function 
of the distance and structural transitions in the ground 
state configurations are examined. 



III. RESULTS 

The results for the ground and metastable states are 
summarised in Table ||. The electron configurations are 
given at four different distances between the artificial 
atoms. The ground state energy and the correspond- 
ing configuration at the four distances is represented in 
the row following the particle number N. If there ex- 
ist metastable states at the given N and d the energy 
difference AE/N to the ground state and the electron 
configuration for the metastable state is also reported. 
However, not all metastable configurations are marked 
in Table ||, since when starting the simulation from ran- 
dom positions more electrons can be trapped in one ar- 
tificial atom than in the other. Only metastable states 
with either the same number of electrons per atom (for 
even N) or only one more at one than the other atom 
(for odd N) are reported. The notation for the configura- 
tions in a single artificial atom is chosen so that electrons 
are thought to be organised in (nearly) concentric shells 
around the potential minimum: (Ni,N2,N3), where Ni 
denotes the number of electrons in the innermost shell, 
N2 the next shell and N3 the number of electrons in the 
outermost shell. (For < 20 only three shells are occu- 
pied). For laterally coupled two-atom artificial molecule 
we have chosen the following notation for configurations: 
At d = 200 A the confi guration is marked as if it would 
still be on a single atom centred around the midpoint 
connecting the two atoms. At c? = 600 and 1000 A the 
configurations are given as configurations of two sepa- 
rated atoms. 

For example, as Table || and Fig ^ (a) indicate, with 
eight particles in the single artificial atom (c? = A) 
the ground state is (1,7), one electron in the centre and 
seven electrons in the circular shell, and there exist no 
metastable states. At c? = 200 A a new ground state 
has appeared with configuration (2,6) (Fig. [l] (b)) while 
(1,7) has changed to a metastable state (see also Table 
|). At distances d = 600 A and d = 1000 A the notation 
is changed to two-atom configurations and for A^ = 8 the 
ground state is marked with (4), (4), see Figs. |l] (c) and 
(e). 

The notation for configurations is not always exhaus- 
tive. The relative orientation of different shells and 
the relative orientations of the configurations of the two 
atoms at d = 600, 1000 A do not always become clear 
from Table ^. For example, when either or both atoms 
are left with four electrons (A = 7, 8, 9), the orientation 
of the square(s) relative to the other atom changes as the 
distance is increased. At smaller distances the position 
of the square of four electrons is such that the tips of 
the squares are in the same line with the positions of the 
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minima (see Fig. I (c) ioT N = 8 at d = 600 A). As 
the distance is increased the square (or two squares with 
= 8) turns onto its side (see Fig. |l| (e) for = 8 
at d = 1000 A). For A^ = 8 at d = 600 A there also 
exists a mctastable state where one of the squares is ly- 
ing on its side and the other on the tip (Fig. |] (d)). 
Even though we divide electrons into shells in our nota- 
tion it does not mean that the shells are strictly circular 
even for d = 0. This can be seen clearly for A^ = 12 in 
Fig 1^ (f), where the outer shell resembles more like an 
incomplete triangle with the tips missing. The configura- 
tion marked with (1,5)* in the two-atom configurations 
in Table | cannot be identified strictly to (6) but neither 
to (1,5). Therefore we choose the notation (1,5)* to de- 
scribe the configuration. The difference between (1,5)* 
and (1,5) can be seen with the (1,5)*, (1,5) metastable 
state in Fig. (g). The configurations of the ground and 
mctastable states for A^ = 13 at d = 200 A are marked 
in the same way in the Table, but are different as can 
be seen in Figs. |^ (h) and (i). For the highest energy 
metastable state for A^ = 17, d — 200 A the two-atom 
notation would have described the configuration better. 
Fig. |l] (1). One metastable state for A^ = 19 and the 
ground state and one metastable state for A^ = 20 at 
d = 200 A could not be described with the shell struc- 
ture notation. The configurations are depicted in Figs. ^ 
(m),(n) and (o), respectively. 

As the distance between the atoms is increased it is 
not always clear whether the electrons just follow the 
two atoms drawn apart and continuously change to two 
separated atoms. Sometimes metastable states change 
to a ground state and the ground state to a metastable 
state as the distance between the atoms is increased. The 
clearest example can be seen in Table | for six electrons 
between d = and d = 200 A. At d = the (1,5) config- 
uration is the ground state and (6) the metastable state. 
At d = 200 A it is the other way around: (6) is the ground 
and (1,5) a metastable state. The energy as a function 
of distance for two alternative configurations is shown in 
Fig.^ (a). The transition point, marked with a small cir- 
cle, is at d = 111.6 A. The transition is continuous with 
respect to energy as a function of distance, but the curva- 
ture of the iJ(d)-curve is different and therefore the first 
derivative of energy with respect to d is discontinuous. 
This is a first-order discontinuous structural transition 
in the electron configuration. Hereafter, by discontinu- 
ous structural transitions we mean the qualitative change 
in the ground state electron configuration which is dis- 
continuous with respect to dE/dd at the transition point. 

In addition to A^ = 6, for A^ = 8, 11, 16 and 19 one 
qualitative change in the electron configuration is ob- 
served as a function of d. For A^ = 8 at d = 135.9 A 
the electron configuration changes from (1,7) to (2,6), 
see Fig |^ (b). Notice that the (2,6) configuration is not 
stable at the limit of one atom and becomes unstable 
approximately below d — 17 A. For A^ = 11 there ex- 
ists one metastable state, (2,9), at d = 200 A, which 
at 214.2 A changes to a ground state as is depicted in 
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(b) 



(2,6) 



N = 8 E/N = 8.850 meV 



(e) 



(4),(4) 



N = 13 E/N = 18.624 meV 




(h) 



(4,9) 



N = 15 AE/N = 0.035 meV 
• • • 

• • • m 



(k) 



• • • 



(5,10) 



N = 20 E/N = 26.843 meV 



(n) 



(-) 



N = 8 E/N = 9.932 jneV 



(C) 



(D,(4) 



N = 12 E/N = 19.111 meV 

(f) (3,9) 



N = 13 AE/N = 0.024 jneV 



(i) (4,9) 



N = 17 AE/N = 0.018 meV 

• • • • 

• •• • • 
• • • • 



(I) 



(6,11) 



N = 20 AE/N = 0.001 meV 
• • • • 



(o) 



(-) 



FIG. 1: (a) - (o). Electron configurations of selected ground 
and metastable states. The energy per particle {E/N) of 
a ground state configuration and the energy difference to 
the ground state per particle (AE/N) for metastable con- 
figuration is given. The configuration is marked on the 
lower right corner. To make it easier to distinguish dif- 
ferent configurations and the distance between the minima 
(d = 0, 200, 600, 1000 A) a circle with 330 A radius is plotted 
around each parabolic potential minimum. A dashed circle 
indicates a metastable state. 



Fig. (c). With A^ = 16 the configuration changes 
from (1,5,10) to (6,10) at d = 102.4 A. The energy dif- 
ferences between ground and metastable states are small 
and metastable states exist only in the proximity of the 
transition point (Fig. | (d)). At d = 80 A the (1,5,10) 
configuration which is plotted below the energy curve is 
still the ground state and the (6,10) configuration plot- 
ted above the energy curve has appeared as a metastable 
state. For A^ = 19 Fig. | (e) shows how the (1,6,12) 
configuration, which is the ground state at d = 200 A, 
has changed to a rather unsymmetric configuration at 
d = 250 A. This unsymmetric configuration changes con- 
tinuously to (2, 8), (1,8) (lower middle plot in (e)) which 
at 586.0 A changes discontinuously to (2, 8), (2, 7) config- 
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TABLE I; Ground and metastable configurations and corresponding energies in meV/particle at four studied distances (d — 
0, 200,600, 1000 A). The energies of metastable states, AE/N, are given as meV/particle above the ground state energy. The 
configurations of electron clusters in the two-atom molecule are described with concentric shells located around the centre of the 
system at d = 0, 200 A and as two separate electron clusters located near the minimum of one of the two atoms at d = 600, 1000 

A. 
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uration. The (2, 8), (2, 7) configuration appears approxi- 
mately at 415 A as a metastable state. The (1,8), (2, 8) 
persists as a metastable state to the greatest studied dis- 
tance of 1000 A as is also indicated in Table |. The other 
metastable states at d = and 200 A do not change to a 
ground state. 

More than one discontinuous transformation in the 
electron configurations is found for = 17 and 18, see 
Figs. I (a) and (b). For iV = 17 the (1,6,10) changes to 
(6,11) a.t d — 145.0 A. The (6,11) configuration parts to 
(1,7), (1,8) two-atom configuration, which is the ground 
state up to d = 501.1 A, where the (2,7) configura- 



tion of the other atom becomes more stable than (1,8), 
thus (1,7),(1,8)->(1,7),(2,7). The (1,7), (2,7) configura- 
tion persists as a ground state to greater distances. Qual- 
itatively similar transformations are seen for = 18 as 
for N = 17. First the centred cluster (1,6,11) changes 
to an open configuration (6,12) at d = 233.9 A. The 
open configuration follows the separation of atoms adopt- 
ing the configuration (1,8), (1,8), where as in = 17, 
the (2,7) becomes more stable than (1,8). However, we 
now see two discontinuous transformations. The first at 
d = 522.8 A when (1,8), (1,8)^(2,7), (1,8) and the second 
at d = 593.4 A when (2,7), (1,8)^(2,7), (2,7). 
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FIG. 2: (a) - (e). Energy per particle as a function of distance for A'' = 6, 8, 11, 16, 19. The small circle indicates the transition 
point. To make it easier to distinguish different configurations and the distance between the minima a circle with 330 A radius 
is plotted around each parabolic minimum. A dashed circle indicates a metastable state. 



For other electron numbers besides the reported N — 
6, 8, 16, 17, 18, 19 we do not observe discontinuous struc- 
tural transitions in the electron configuration as the dis- 
tance between the two atoms is increased. A few ex- 
amples of continuous electron configuration changes are 
shown in Fig. ^. For = 12 the (3,9) configura- 
tion transforms continuously to resemble the (6), (1,5)* 



two-atom configuration. Between 200 A and 450 A the 
row of electrons pushes itself forward when the atoms 
move apart, resulting in the symmetric configuration 
(1,5)*, (1,5)*. For TV = 14 the electron configuration 
follows the separation of atoms in a symmetrical form, 
but after d — 250 A both sides start to twist towards 
the (1,6), (1,6) two-atom configuration. For = 20 the 
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FIG. 3: Energy as a function of distance for TV = 17 and 
A'' = 18 with the ground state configurations along the energy 
curve. The small circles indicate the discontinuous structural 
transition points. 



transformation is hard to describe, but it is continuous. 
The three distinct metastable states at d = and 200 
A, seen in Table | never change to a ground state and 
vanish at other distances. 

The changes in energy per particle of the ground states 
at the four studied distances as a function of N are shown 
m Fig. |. At d = there are small troughs at iV = 3, 6, 10 
and 17, at adding the fourth, seventh, eleventh and eigh- 
teenth particle. Moving to greater distances between the 
atoms, the change in the chemical potential is clearly 
peaked. Going to an odd number of particles increases 
the chemical potential much more than going to an even 
number of particles. Interesting is the intermediate dis- 
tance oi d — 200 A where this trend is observed for 
N — 2,3,4 and 9, 10, 11, but otherwise the curve shows 
no clear structure and does not strictly follow the shape 
of the d = A curve either. 




FIG. 4: (a) - (c) Ground state electron configurations along 
E{d) curve for iV = 12, 14 and 20. 



IV. DISCUSSION 

At d = (single artificial atom) our results are in 
agreement with other Monte Carlo (MC) and molec- 
ular dynamics (MD) studies with ■ uaijabpUc ■ ^anfme- 
ment and pure Coulomb interactip«.tilE3E3EJEj'E3 How- 
ever, for = 17 Bolton et al. lH3 obtain the (1,5,11) 
configuration instead of (1,6,10) which our and other 
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FIG. 5: Change in the chemical potential {En+i/{N + 1) — 
En/N) at d = 0, 200, 600 and 1000 A. 

calculationsilii0'0 predict. Ref. ^^-jnay contain an 
error since in a later work by BoltonEil the configura- 
tion for = 17 was reported to be (1,6,10). There is 
also a difljescnce for = 21 in Ref. ^ which was later 
corrected.Eil 

Besides calculating the ground state configurations 
Kong et al. EiJ also examined metastable states for N = 
1—40. Our results are in agreement (we calculated config- 
urations only for N < 20) both in ground and metastable 
states except that for = 18 Kong et al. found two 
metastable states whereas we see only one. In addition 
to the (1,6,11) and (1,7,10) they also obtain (6,12) as a 
metastable state. We repeated the simulation with 3000 
independent test runs, but were still unable to find the 
(6,12) configuration. 

We can conclude that different calculations for the 
confinement and l/r interaction potential are in good 
general agreement. The few experiments on charged par- 
ticles trapped in 2D as well as calculations with different 
forms of interaction and confinement potentials reveal 
also different configurations for the cluster patterns. The 
interaction between the particles could be logarithmic, 
which is the case with infinite charged lines moving in 
2D (vortex lines etc.) or perhaps Yukawa type with a 
strong but short-range repulsion (screened Coulomb in- 
teraction). The form of the .confinement is usually chosen 
to be parabolic (Lai and IcJ tested also a steeper con- 
finement with r'* contribution). However, for the ques- 
tion whether the potential in experiments with clusters is 
parabolic there is no clear answer. Therefore it is not sur- 
prising that the experiments and also calculations with 
different functional forms of interaction and confinement 
result in different cluster patterns. 



Lai and calculated and summarised the configura- 
tion patterns with different interactions and tested also 
contribution to the confinemept and compared thejjesults 
to dust particle experiments. l3 Saint Jean et al. E3 mea- 
sured the configurations with electrostatically interacting 
charged balls of millimeter size moving on a plane con- 
ductor. They made a comparison with simulations and 
quite surprisingly found the best agreement with a rela- 
tively old simulation with vortex lines in a superfiuiclc£l 
with logarithmic interaction, which again mas not in 
agreement with the dust particle experimentstS nor with, 
the purely logarithmic interaction used by Lai and I.u 
Despite the differences there are some particle num- 
bers where the configuration seems to be the same re- 
gardless of the experiment or functional form of the in- 
teraction or confinement. These particle numbers are 
N = 3,4,5,7,10,12,14,19. 

Partoens et alS3 examined the ground states of even 
number of classical electrons evenly distributed in two 
vertically coupled artificial atoms as a function of the 
distance between the atoms. As in our study of laterally 
coupled atoms discontinuous transitions between config- 
urations occur as a function of the distance d. The dif- 
ference is that for vertically coupled atoms one can see 
intuitively that transitions should occur between d = 
and c? — > oo whereas in laterally coupled atoms the con- 
figurations can be pulled apart with some N without 
qualitative (discontinuous) changes in the electron con- 
figurations. For vertically coupled atoms discontinuous 
transitions (first order with respect to energy) in electron 
configurations are observed for all even < 20 whereas 
for laterally coupled atoms we see also purely continuous 
changes as c? — > oo. 

To summarise, we have calculated ground and 
metastable state configurations of classical point charges 
confined in two dimensions with two laterally coupled 
parabolic potential wells. Ground and metastable elec- 
tron configurations were studied as a function of the dis- 
tance between the atoms and discontinuous (in dE/dd) 
transitions in the ground state configurations were ob- 
served for particle numbers A^ — 6,8,11,16,17,18,19. 
The configurations of purely classical electrons in lat- 
erally coupled two-minima potential have an interesting 
and complex spectrum as the distance between the min- 
ima is changed. 
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